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We study the integrability of gravity-matter systems in D — 2 spatial dimensions 
with matter related to a symmetric space G/K using the well-known linear systems 
of Belinski-Zakharov (BZ) and Breitenlohner-Maison (BM). The linear system of BM 
makes the group structure of the Geroch group manifest and we analyse the relation 
of this group structure to the inverse scattering method of the BZ approach in general. 
Concrete solution generating methods are exhibited in the BM approach in the so- 
called soliton transformation sector where the analysis becomes purely algebraic. As a 
novel example we construct the Kerr-NUT solution by solving the appropriate purely 
algebraic Riemann-Hilbert problem in the BM approach. 
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1 Introduction 

In late 1970s and early 1980s a large variety of solution generating techniques for the four- 
dimensional vacuum Einstein equations and Einstein-Maxwell equations were explored, for an 
overview see [1]. It was later realized in the case of two commuting and hypersurface orthogonal 
Killing vectors that all these approaches are nothing but different manifestations of the integra- 
bility of the corresponding effectively two-dimensional system that can be exhibited by means of 
a linear system or Lax pair [2, 3, 4]. Several authors made efforts to find interrelations between 
these different methods. Cosgrove [5] took a computational approach whereas Breitenlohner and 
Maison [3] concentrated on unraveling the group theoretical structure behind these techniques, 
taking up ideas of Geroch [6] and Julia [7]. Other relevant work includes [8, 9, 10, 11, 12, 13, 14]. 
In the present work, we revisit these studies to further elucidate the interrelation between the 
various techniques. 
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The best developed technique for generating solutions is that of Belinski and Zakharov [2, 4], 
henceforth BZ for short and often called the inverse scattering method. It has been very suc- 
cessful in constructing novel solutions in both four- and five-dimensional vacuum gravity. The 
method involves some rather special adjustments to certain quantities before new physical so- 
lutions can be obtained. However, there are problems when applying this method to different 
gravity-matter systems like those of interest in supergravity where the implementation of the 
same adjustments fails and it is not guaranteed that an inverse scattering transformation pre- 
serves all features required of a solution to the gravity-matter equations [15]. In the group 
theoretical framework of Breitenlohner and Maison (BM) [3] this problem does not arise since 
the solution generating transformations form the so-called Geroch group (an affine group) and 
by the group property any transformation will generate a new solution. The drawback of the 
BM method is that it is not as easy to implement and does not always operate directly on the 
physical quantities. Despite these shortcomings, the promise of the BM approach is that it can 
be applied to various general settings of interest. 

In order to illustrate this point, we consider for concreteness D = 4 gravity with a space- 
like and a time-like Killing vector such that the system is effectively two-dimensional and we 
are in the realm of so-called stationary axisymmetric solutions. The infinite-dimensional affine 
symmetry group can be viewed as the closure of two finite-dimensional symmetry groups that 
act on the space of solutions [6, 3]. The first one is typically called the Matzner-Misner group 
SL(2, 1)mm and consists of area preserving diffeomorphisms of the orbit space of the two Killing 
vectors (the 'torus' that one reduced on, the Killing vector orbits are not necessarily compact). 
The other group is called the Ehlers group SL(2,R)e and is a hidden symmetry already of the 
three-dimensional model. The fields that it acts on are partly formed from dualising a Kaluza- 
Klein vector field (in D = 3) to a scalar field; therefore it does not act directly on the metric 
components. Combining SL(2,M)mm and SL(2,M)e yields the infinite-dimensional affine Geroch 
group [6, 3]. 

The linear systems of BZ and BM display the integrability of the D = 2 theory in terms of 
generating functions that depend on spectral parameters besides the dependence on the two- 
dimensional space(-time) coordinates x m . There are in principle four choices for the linear 
system and generating functions at this point that can be summarized in the following diagram: 





Ehlers 


Matzner-Misner 


BZ 


*e(A,x) 


< — > *mm(A,x) 




t 


t 


BM 




< — ► Vmm(*,z) 



The horizontal arrows in this diagram correspond to duality transformations between some of 
the two-dimensional scalar fields. At the level of the Lie group this corresponds to the action 
of an (outer) automorphism, now known as the Kramer-Neugebauer mapping [16, 3, 17]. The 
vertical lines relate the two linear systems and are the main subject of the present work. Note 
that the spectral parameters of the generating functions are different in the two rows. The 
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traditional inverse scattering method of BZ operates on the physical variables and so is located 
in the Matzner-Misner column of the diagram. In fact it uses a slight extension of the Matzner- 
Misner variables in that it also treats the volume of the two-dimensional orbit space and so 
corresponds to GL(2,IR)mm rather than SL(2,M)mm- In our work, we will focus mainly on the 
left column where the group theoretical structures are a bit nicer, e.g., empty Minkowski space 
corresponds to the identity function. We will always display the subscripts E and MM to avoid 
any confusion. 

It follows from the work of BM that any gravity-coupled matter system in D = 3 with 
symmetric space target Ge/Ke becomes integrable when an additional Killing symmetry is 
present. Here, Ge is the generalization of the Ehlers symmetry SL(2, R)e to other systems. The 
group of symmetry transformations form the (non-twisted) affine extension of Ge and we will 
restrict to simple Ge for simplicity. Beyond this general fact, it is quite hard to construct actual 
solutions explicitly using this method since one has to solve a matrix valued Riemann-Hilbert 
problem [1, 3, 18], sometimes also referred to as Birkhoff factorization. In the case analogous to 
the soliton transformations of BZ, however, the problem reduces to a linear algebra system that 
can be solved as shown in [19] which we will review below since the work is not published. 

Our chief motivation to study the different formulations of integrability is to find their 
interrelation and, by this, to make new methods available for generating solutions of (super- 
gravity beyond the cases that have been covered so far. Of particular interest are cases when 
Ge is an exceptional symmetry group and other cases that arise in string theory. There are 
no known established systematic techniques for constructing non-supersymmetric solutions that 
exploit the integrability structure of supergravity theories. With such techniques at hand one 
can construct a variety of new solutions, in particular, new black hole solutions generalizing 
[20, 21] and new fuzzball solutions generalizing [22, 23]. 

The structure of this article is as follows. In sections 2 and 3 we introduce the linear 
systems of Breitenlohner-Maison and of Belinski-Zakharov in Ehlers form and elucidate their 
interrelation. In section 4 we study meromorphic generating functions and solve the related 
linear systems algebraically. Section 5 gives the explicit example of the Kerr-NUT solution in 
both cases. We conclude in section 6. Certain technical computations have been relegated to 
the appendices. 

2 Breitenlohner-Maison linear system 

In this section we review the linear system of Breitenlohner and Maison (BM) [3, 24, 25, 26]. 
This linear system arises from considering a D = 3 gravity-matter system with target being 
a symmetric space Ge/Ke where Ge is the global symmetry group that we refer to as Ehlers 
symmetry. Our conventions are such that the three-dimensional space has signature + + + 
and we consider the space of solutions admitting an additional (axial) symmetry so that the 
theory becomes effectively two-dimensional. This two-dimensional system can be shown to be 
integrable. 
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2.1 D = 3 model with Ge/Ke matter 

Consider a three-dimensional Euclidean gravity-matter system with action 

Se = J d 3 x^ (r - \{P^\Pv,u)g^ , (2.1) 

where (-|-) is the symmetric invariant bilinear form on the Lie algebra gE of the real Lie group 
Ge- The group Ke is a subgroup of Ge with dimension equal to that of the maximal compact 
subgroup and the coset space Ge/Ke is a (pseudo-)Riemannian symmetric space. The subgroup 
Ke is fixed by some involutive automorphism r on Ge [27]. This induces an involution 9 on 
the Lie algebra (equal to the Cartan involution when Ke is maximally compact) which splits 
0E = %e ffi Pe- Let Ve £ Ge/Ke be a coset representative (e.g. in Borel gauge according to the 
Iwasawa theorem in a patch where it applies); then we decompose 

OvVeVe 1 = Pe,„ + Qe,h (2.2) 

with 

Qe,[m = \ (Q^VeVe 1 - id^VEV^)*) , (2.3a) 

Pe^ = \ (d^E 1 + (9„VeV e x )#) , (2.3b) 

where we have defined the 'generalized transpose' = —6(x) on the Lie algebra. (For s[(n,IR) 
it is the standard transpose when t = so(n).) The two components in the splitting (2.3) satisfy 

Qf = ~Qe, Pt = Pe- (2.4) 

The map # is an anti-involution and its group version will be denoted by the same symbol and 
is also an anti-involution: For g,h 6 Ge one has g& = r(g^ 1 ) = r(g)~ l and (gh)& = bfig& . 
The symmetry transformations acting on Ve are 

Ve(x) -> kE(x)V E (x)g E (2.5) 

with constant gE £ Ge (global transformations) and varying &e £ Ke (gauge transformations). 
Under these transformations, the Lie algebra valued quantities transform as 

Qe,ii -> kEQE^k^ 1 + d^kEk^ 1 , (2.6a) 
Pe,» -> k E QE,^\ (2.6b) 

i.e., Qe as a KE-connection and Pe transforms KE-covariantly. 

A useful quantity associated with Ve £ Ge/Ke is the 'monodromy matrix' 

M E = V*Ve (2.7) 

that transforms as Me — > g^MEgE under (2.5) since k^kE = 1. It is hence insensitive to 
the KE-gauge chosen for the coset representative Ve and only transforms under the global Ge 
transformation. 
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The D = 3 equations of motion derived from (2.1) are 

R»u ~ ^(PeJPe,v) = 0, (2.8a) 

d„ {V=9g^Pk,u) - V=gg^ [Qe,„, Pe,„] = o. (2.8b) 

For convenience, we introduce the i^E-covariant derivative 

£>E,„ = d„- [Qe, m , •] (2.9) 



in terms of which (2.8b) becomes 



DeAV^Pe) = 0- (2-10) 



2.2 Reduction to D = 2 

If the system admits an axial isometry d$ we reduce the metric according to 

dsl = fldsl+p 2 d<p 2 . (2.11) 

The function /e will be referred to as the conformal factor of the effective two-dimensional 
metric. We label the two-dimensional coordinates as x m . The two-dimensional metric ds 2 . is 
assumed to be flat by appropriate choice of coordinates. Note that there is no Kaluza-Klein 
vector A m of the type [d(j> + A m dx m ) 2 in (2.11) since it carries no degrees of freedom and can 
be set to zero without loss of generality. 

With this ansatz, the equations of motion (2.8) then imply 

□p = 0, (2.12) 

which we solve by choosing Weyl canonical coordinates x m = (p, z) on the flat two-dimensional 
space so that ds 2 , = dp 2 +dz 2 . We let z be the conjugate variable to p such that d p = e pz d z = *2<9 2 . 
It is often useful to combine the two real variables into a single complex variable (and its complex 
conjugate) which we denote by 

x ± = ^(z T ip) (2.13) 

in analogy with light-cone coordinates that would arise if the two-dimensional base was Minkowskian 
as for colliding plane wave solutions. Note that (2.13) implies d±p = ±i and *2<9± = ±id±. 
With this choice the remaining equations are equivalent to 



D m (pP£) = 0. (2.14b) 



±if^d±fE = ^(Pe,±\Pe,±), (2.14a) 



The first equation is a constraint and yields the conformal factor by a single integration and is 
therefore of secondary interest. The main equation of interest in this paper is the last equation 
(2.14b). A key result of [3] is that this equation is integrable and has an underlying symmetry 
structure associated with the affine extension of Ge- 
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2.3 BM linear system in general 



Consider the generalization of Ve £ Ge/Ke to also depend on some spectral parameter t (we 
suppress the space-time dependence) V E — > Ve(£), where 



^ 0) +< ) + ^ + 



importantly is regular in t around t = and the limit t — > gives back the original Ve- 



limVE(t) = Ve(0) = V£ )) :=V e . 



Consider the linear system [3, 26] 



W 1 = ^% + Oe,± ) 
1 i it 



(2.15) 



(2.16) 



(2.17) 



where Qe and -Ffe are independent of the spectral parameter t and are defined in terms of the 
t-independent Ve as in (2.3). 1 The integrability condition for (2.17) is equivalent to the equation 
of motion (2.14b) if and only if the spectral parameter t satisfies the differential equation 



t L d±t = -——.P l d±p. 
1 ± it 



(2.18) 



This differential equation can be integrated [3] to an equation for t which we write in the more 
conventional form in terms of the Weyl coordinates (p,z) (cf. (2.13)) 



2 2t , 

r (z-w) -i = o, 

p 

where w is an integration constant. This quadratic equation has two solution branches 

1 



t± 



P L 



(z — w) =b \J (z — w) 2 + p 2 



(2.19) 



(2.20) 



Equation (2.19) defines a two-sheeted Riemann surface over the two-dimensional flat base. We 
take the solution with the + sign to be the physical sheet and when we write t we always mean 
t + unless indicated otherwise. We will refer to t as the space-time dependent spectral parameter 
and to w as the constant spectral parameter. 2 

We will refer to (2.17) as the BM linear system (in Ehlers form) and we have just reviewed 
how its integrability condition gives rise to the equations of motion of the D = 3 gravity-matter 



1 In a general coordinate system, the linear system takes the form 



PE.m - 



2t 



l + t 2 "" l + t 2 

2 When written in terms of the 'light-cone' coordinates (2.13), equation (2.20) becomes 

t± = 



s/w — 2x + =p \Au — 2x 



\/w — 2x+ ± \Au — 2x 
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system (2.1) in the presence of a Killing isometry. This establishes the integrability of the 
equation (2.14b); the conformal factor /e can then be obtained by integrating (2.14a) [3]. The 
linear system (2.17) is vastly underdetermined since it represents two differential equations for 
a function of three variables. There is an infinity of integration constants associated with this 
system. 

Besides giving the integrability of (2.14b) the BM linear system also serves to unveil the 
group theory underlying the system. The original t- independent Ve G Ge/Ke transformed 
under global Ehlers transformations gE G Ge as Ve — > &eVe<7e ( c f- (2-5)), where &e G Ke 
is the usual local compensator required to restore a chosen gauge for the coset representative. 
The presence of the spectral parameter now suggests to enlarge the set of global symmetry 
transformations by allowing qe to depend on the constant spectral parameter w: 

Ve(*) ->• k E (t)VE(t)gE(w). (2.21) 

As indicated, the compensator is now also i-dependent as it has to be chosen such that the 
transformed Ve(£) is regular around t = as in (2.15). This enlarged set of global transformations 
consists therefore of functions #e(^) ; he., maps of the type C — > Ge, where we impose that gE(w) 
admits an expansion around w = oo in order to remain expandable as in (2.15). These maps 
include transformations from S 1 C C into Ge and (under additional regularity assumptions) 
will lead to the loop group Ge associated with the Ehlers group Ge- Therefore the group 
underlying the integrability in D = 2 includes the infinite-dimensional loop group; in fact the 
extension to the full afline group is active [7, 3] where the central extension acts on the conformal 
factor /e (see below). We note that besides the afline group one can also define the action of 
the (centerless) Virasoro algebra which arises from arbitrary reparametrisations of the constant 
spectral parameter [28, 29, 24]. Together with the infinitesimal afline transformations one obtains 
a semi-direct product in the standard way. We will not use the Virasoro symmetry in this paper. 
The involution # extends to functions Ve(£) by 

(VeW) # = V* (-i) . (2.22) 

One can use # to split the loop algebra into an invariant and an anti-invariant part (generalizing 
pE and €e above). Now, it is important that the right hand side of (2.17) is anti-invariant under 
(the Lie algebra version of) # and therefore belongs to the 'compact' subalgebra of the affine 
algebra based on qe- This anti-invariance implies that ((Ve^))*)" 1 is a solution to the linear 
system if Ve(£) is a solution. In general, the two solutions related by this involutive mapping 
will be different. The mapping implies that the monodromy matrix 

Me = M E (w) = (V E (t)) # V E (t) = V* VB(t) (2.23) 

is independent of the space-time coordinates x m and therefore a function of w alone. The matrix 
A^e(^) is invariant under the application of # and simultaneously exchanging t — > —1/t. We 
note that this is evident from (2.19) which implies that w(t,x) is invariant under t — > —1/t. 
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Constructing new solutions of the linear system by means of the Geroch symmetry proceeds 
along the following chain of steps 3 

V E V E (t) MeH -> M 9 E (w) ->• V|(t) ->• y E 9 , (2.24) 

where we introduced the notation 

.MeH -> >*eW : = 9*(w)MeH9e{w) (2-25) 

for the transformed solution. The individual steps in (2.24) starting from a given seed solution 
Ve are: (i) find a corresponding generating function Ve(£) that solves the linear system (2.17), 
(ii) compute the associated monodromy (iii) transform the monodromy under a global trans- 
formation g(w) as in (2.25), (iv) factorize the new monodromy into a new generating function 
Vfi(t) and (v) take the limit t — > to find the new solution. 

For practical purposes, the main difficulty resides in step (iv) in factorizing the transformed 
M%(w) as 

Mi(w) = (Vg(t))#Vg(t) (2.26) 

with the new V^(t) having an expansion as in (2.15). This is a Riemann-Hilbert problem [3] 
whose solution is in general hard to obtain. In the particular case of meromorphic M^(w) with 
single poles in w of certain simple type one can reduce the problem to a set of linear algebraic 
equations. This is the case of soliton charging transformations that will be discussed further in 
section 4. Once the new V^(t) has been obtained, one can recover the solution to the gravity- 
matter system (2.14) by taking the limit i — ?• and obtain VS 6 Ge/Ke that characterises the 
physical fields. 

Besides the knowledge of the coset 'scalars' V E E Ge/Ke one also requires the new conformal 
factor /| in (2.14). This can be obtained from a simple integration of (2.14a) but it also follows 
from group theoretic properties using the central extension. This is discussed in detail in [3] to 
which we refer for the general expression. In section 4 we will present the formula in the case of 
soliton transformations. 

We note that a trivial solution of the equations (2.17) and (2.14a) is given by 

Ve(*) = 1 and / E = l. (2.27) 

This solution will be referred to as flat space as it corresponds to the Minkowski vacuum in the 
four-dimensional case. 

3 Belinski— Zakharov linear system 

In this section we present the linear system used by Belinski and Zakharov (BZ) [2, 4]. We will 
not present it in the standard form which uses what was called the Matzner-Misner formulation 

3 We present here the solution generating method based on the monodromy matrix A4e(w). Alternatively, 
one could work at the level of the generating function Ve(£) and the transformation (2.21); however, the step of 
finding the compensator k^it) in (2.21) is typically very hard. 
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in the introduction. Rather, we will use the Ehlers description to make contact with the discus- 
sion in the preceding section. (The Matzner— Misner version and its relation to BM is discussed 
in appendix B.) 



3.1 BZ Ehlers linear system 

Equation (2.14b) for the Ge/Ke coset fields admits an alternative Lax pair that can be written 
as 

X z + p z X z + p z 

where A is the (space-time independent) spectral parameter of BZ and ^e(X,p, z) is the gener- 
ating function such that the matrix Me = Ve of (2.7) is recovered for A = 0: 

M E {p,z) = y E (0,P,z). (3.2) 

The matrices U,V are defined as U = pd p M E M^ , V = pdzM^M^ 1 , and the differential 
operators D\,D<i are 

X z + p z X z + p z 

The operators -Di and Z?2 commute and the associated integrability condition of the linear 
system (3.1) is equivalent to the desired non-linear equation (2.14b). 

Solutions of the BZ linear system (3.1) are constructed using the inverse scattering method [2]. 
One starts from a 'seed' VPe.Oj that is 'dressed' to obtain a new solution ^e through 

*e(A) = x(A)*e,o(A) (3.4) 

where x 1S called the dressing matrix and it depends on the spectral parameter A. The seed ^e.o 
corresponds to a solution of (3.1) for a seed solution Me,o- We can take it to be the identity 
matrix ^e,o = 1> which corresponds to taking the seed solution to be flat space 4 . In order for 
the 'dressed' \I/e to also solve the linear system (3.1) the dressing matrix has to satisfy its own 
linear system 

pV-XU pV - XU pU + XV pU + XV 

DlX = ^Tp^ X ~ X A2 + p2 ' D2X = l^Tp^ X ~ X A2 + P 2 • ( 3 - 5 ) 

In addition, the matrix x must satisfy further constraints in order to ensure that the new 
solution Me{p,z) = ^(0,p,z) is real, satisfies M^ = Me and is a representative of the coset 
Ge/Ke [2, 15, 30]. 



4 Note that this differs from the more common (and compficated) expression for fiat space in the Matzner- 
Misner form [2, 4]. 
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3.2 Relation between the two linear systems 

Compared to the discussion of the BM linear system, the differential operators D\ and Z?2 of 
(3.3) can be demystified by thinking of A as space-time dependent, so that [3] 

Di = d z = 9 Z I a fixed + d z \d x , D 2 = d p = d p \xR xc d + dp\d\. (3.6) 

If the spacetime dependence of the spectral parameter A is given by 



one recovers (3.3). The solution to (3.7) is given by 



A(p,z) = [w-z)t + (3-8) 

where w is an integration constant. Comparing to (2.20), it follows that from this viewpoint the 
relation of the BZ spectral parameter A to t in the BM approach [3, 15] is 5 

\(p,z) = -pt(p,z). (3.9) 

The relation between the two generating functions Ve of (2.17) and ^e of (3.1) is given by 

* E (A,x) = V*(x)V E (t,x) (3.10) 

where one also has to use (3.9). Note that on the right-hand side we have once the spectral 
parameter independent Ve(x) = Ve(0, x) and once the full VE(t,x). This obscures the action 
(2.21) of the full affine Geroch group since the transformation of Ve(x) under affine elements is 
complicated. 

In the following we restrict to Ge = SL(n,R) for concreteness. In that case Me has to be a 
symmetric matrix. For other groups, there will be different conditions on some of the quantities 
introduced below. 

If the matrices Me and Me,o obtained by the A — > limits of ^e and ^e.o are symmetric, 
then 

-,2 - 



X / (A)=Mex T_1 (-j)m^ , (3.11) 



solves exactly the same linear system (3.5) as x [2]- Given this observation, one has that x'(^) 
is related to x(A) through some arbitrary matrix C(w) via 

x'(A)* E ,o = x(A)*e,oCH. (3.12) 

This reflects the fact that the linear system (3.5) is underdetermined and C(w) corresponds to 
a gauge freedom of (3.1). However, Belinski and Zakharov demand 

x'(A)=x(A), (3.13) 



3 The sign in this relation in [3] appears incorrect. 
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which corresponds to fixing the gauge freedom of (3.12). In addition, they do not require x(A) 
to satisfy the coset constraint detx(A) = 1. Since detx 7^ 1 one has that the new matrix Me 
does not have unit determinant and so does not represent a physical configuration. To obtain 



Me appropriately. 

Condition (3.13) automatically ensures that the final Me is symmetric, but it is a rather 
strong assumption. Relation (3.13) is a central equation in the BZ inverse scattering framework 
and it fixes an infinite ambiguity in the dressing matrix that corresponds roughly to the Borel 

part of the Geroch group. In other words, demanding (3.13) in the BM framework amounts to 

(n) 

choosing finely tuned integration constants for all the dual potentials Vg with n > 1 introduced 
through (2.15). The transformation (3.11) is similar to the one discussed below (2.22) in the 
BM framework. There one normally does not fix this freedom and so there is no direct analogue 
of (3.13) in the BM approach. 

Due to the complications of detx 7^ 1 and the issues mentioned around (3.10), it is imprac- 
tical to find a satisfactory embedding of the full BZ solution generating technique in the Geroch 
group 6 . The best one can do is to find a representative relation between the BZ-generating 
function ^/e(A) and the group-theoretic BM generating function Ve(£)- This relation is precisely 
equation (3.10) for the Ehlers coset and is obtained in Appendix B for the Matzner-Misner 
coset, see (B.13). 

3.3 Solitonic solutions 

So-called solitonic solutions for the matrix Me correspond to a dressing matrix x(A) with simple 
poles in the complex A-plane. The general iV-soliton solution is obtained by dressing the seed 
solution with a matrix x °f the form 



The matrices Rk and the pole trajectories are functions of p, z only. For each soliton, there 
exist two possible solutions for the pole trajectory /ij. 



where the parameters Wk may generally be complex but for the examples considered here we 
will take them to be real. The pole trajectories with a "+" sign are referred to as solitons and 
the ones with a "— " sign as antisolitons. 

6 Naively one might conclude from (3.12)-(3.13) that it simply corresponds to taking C(w) = 1 from the 
Geroch group point of view. However, this interpretation is not correct. This is because the dressed BZ matrix 
'J'e(A) = x(A)*Pe,o does not directly give the physical matrix M^ hys \ In order to have an interpretation of C(w) 
in the Geroch group, one first needs to construct T^ phys > (A). A procedure to do this was suggested in [5]. Requiring 
something like x (phys) '(A) = X Cphys) (A) will indeed be more amenable to the group theoretic interpretation, but it 
is not the BZ technique. 



the 'physical' matrix M E that fullfills the determinant condition, 



one rescales the matrix 




(3.14) 




(3.15) 
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In order to construct the iV-soliton dressing matrix, one needs to parametrise the residue 
matrices Rk- Here, one has the freedom of introducing certain arbitrary constant parameters 
m ^ob ( w ith b = 1, . . . , n when ^ is represented as an n x n-matrix) for each soliton pk as follows 7 . 
Defining 

and the symmetric matrix Tbz as 



tfiU^.P.*) , , (3-16) 



(r \ m { a k \M Efi ) ab mP 

(fBZjfcZ = 5— , O 1 ') 

P + PkPl 

the elements of the residue matrices Rk are given by 

... (rzk),.7r$ (Meo) 1, 
(Rk) ab = mW £ 1 BZhk K fi ' cb . (3.18) 

The new matrix M^(p, z) = x(0, p, z)\I/e,o(0; p, 2) now reads 

(M E ) afe - (M E , ) ab - ^ — • (3.19) 

The symmetry of this expression is ensured by (3.13). A problem that arises at this stage is 
that possibly the new matrix M E does not satisfy the coset constraint detMs = 1, i.e. is not 
an element of the group SL(n,M)E- In fact the determinant of the new matrix is given by 

detM E = {-l) N p 2N (jl^ detM E,o- (3.20) 

In order to obtain an iV-soliton solution that remains in the group SL(n, K)e, the new matrix 
Me must be multiplied by an overall factor 8 

The overall sign in this expression should be chosen in order to ensure the right metric signature. 
Thus obtained M E phys ' fulfils the constraint detM E phys ^ = 1. Finally, following the discussion in 
[2], the conformal factor for the dressed solution can also be obtained. We find for SL(2,R)e 

N 



(/E phys) ) 2 =^z-^-(n 



N \ N 



Pk 
k=l J 



{Pk ~ Pi) 



detr BZ • flo, (3-22) 



Pk - Pi" 2 

k,l=l, k>l 

where /cbz is an arbitrary numerical constant. For SL(n, K)e a similar but more complicated 
expression can also be written [4]. However, note that for n > 2 the rescaling (3.21) introduces 
fractional powers of p from (3.20) that typically lead to singular solutions. For this reason it is 
more useful to employ the so-called Pomeransky trick [31] for n > 2. In this approach one can 
write a general expression for the conformal factor valid for n > 2 [31, 15]. 



7 The normalization of each of the vectors is arbitrary. Rescaling them by arbitrary constants does not 

change any of the final expressions. 

8 This formula differs slightly from the standard expression in [2] since we are working in the Ehlers description. 
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4 BM Soliton transformations 



In this section we present an algebraic method of generating new solutions of the BM linear 
system (2.17) from a given seed solution. Our discussion closely follows that of [19], see also [32]. 

The method makes use of the constant group element <?e(V) of the Geroch group. We take the 
seed solution to be flat space (2.27) since it is believed that the Geroch group action is transitive 
on the space of solutions and all solutions are related to flat space [18, 3, 13]. As mentioned 
around (2.26), the action of the Geroch group generally leads to a matrix valued Riemann- 
Hilbert problem. In the case when the matrix functions to be factorized are meromorphic in the 
spectral parameter w, the problem can be solved algebraically. This is the case that we focus 
on and we term it the solitonic case. There are a number of (formal) similarities and at the 
same time a number of differences (in details) with the procedure of Belinski and Zakharov [2] 
that we briefly reviewed in section 3.3. In this section we also restrict ourselves to the Ehlers 
SL(n, R)e of D = 2 + n vacuum gravity. In this case the generalized transpose # (at the group 
level) simply becomes the usual matrix transpose. 

4.1 Riemann Hilbert problem 

In section 2 we presented the construction of Me(w) starting with Ve(x). We start with Ve(x), 
solve the linear system (2.17) to find Ve(£, x), and then construct Me(w). We now ask, following 
the steps of (2.24), if we can reverse the process and reconstruct Ve(x) from Me(w), i.e., we 
solve the Riemann-Hilbert problem to factorize Me{ w ) as in (2.23). It is not guaranteed that 
for a general Me{w) such a factorization exists. In fact, one can construct explicit examples 
where it does not exist. This is however a technical problem that will not concern us here. 
Certain aspects of this have been studied in the literature, see e.g. [18]. We only work with 
those matrices Me{w) for which the Riemann-Hilbert problem admits a solution. 

Let us start with a real symmetric unit determinant matrix A4e(w) assuming suitable ana- 
lyticity properties. In particular we assume .Me(oo) = 1. We wish to factorize it as 

Me(w) = A T _{t,x)ME{x)A + {t,x), (4.1) 

with A-(t,x) = A + (— \,x) and Me{x) symmetric and real. Moreover, we require 

det A±(t,x) = 1. (4.2) 

Next, we factorize Me(x) as Me(x) = (x)Ve(x) with a triangular matrix Ve(x) to obtain 

V E (t,x) = V E (x)A + (t,x). (4.3) 

The factorization problem (4.1) can be viewed in two different ways, (i) We start with an 
appropriate Me(w) and solve for Ve(£, x), (ii) we start with a seed VE(t,x) and act with an 
element <te(u>) and attempt to determine the transformed V E (t,x). We take the first viewpoint 
in what follows. 9 It is, however, convenient to have the second viewpoint in the back of one's 

9 The second viewpoint was taken in [19], however, some of their assumptions about pole structures seem too 
restrictive to make that method directly applicable to interesting solutions. 
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mind and relate it to the first one by taking the seed to be flat space. At the level of equations 
this means 



V 9 (t,x) = Z 9 (t,x)V E (t,x) = Z 9 (t,x) -l = Z 9 (t,: 



(4.4) 



where Z+(i,x) is a triangular 'dressing' matrix in the sense of (2.15) that is determined by the 
Geroch group element g E (w). In terms of the monodromy matrix one similarly has 



M 9 E (w) = Vl 



t'' 



Z 9 



Z 9 (t,x)V E (t,x). 



(4.5) 



Introducing 



Z9(t,x) 



Z 9 



Z 9 (t,x), 



one also has 



Z 9 (t,x) 



1 



MUw) [V E (t,x)]- 1 



l-Ml{w)-t = M 9 {w). 



(4.6) 

(4.7) 
(4.8) 



When V E (t,x) ^ II one should take Z 9 (t,x) in (4.7) to be the left hand side of equation (4.1) 
and solve the corresponding factorization problem. In this paper we always work with flat space 
(2.27) as seed solution. Consequently, for notational simplicity we drop the superscript g from 
now on and just think of being given a monodromy M E (w) that needs to be factorized as 
in (4.1). 

4.2 Multisoliton solutions 

The factorization problem can be solved algebraically when the matrix functions to be factorized 
are meromorphic. We now present this factorization explicitly by adapting [19]. We assume that 
M. E (w) has simple poles with residues of rank one. Since det Ai E (w) = 1 and A4e(oo) = 1 the 
inverse matrix A4^ 1 (w;) also has poles at the same points with residues of rank one. If we have 
N poles at points Wk with k = 1, 2, . . . , iV we can express A^e(^) and A / Jg 1 (tf) in the form 



v 



.MeW = a + 

k=l 
N 



Ai 



k=l 



W -w k 

B k 

w-w k ' 



(4.9a) 
(4.9b) 



with symmetric (since M E is symmetric) and constant residue matrices Ak and Bk of rank one. 
This means that we can factorize these matrices as the outer product of vectors 



A k = a k a k a k 



hPkb 



T 
k ' 



(4.10) 
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One could absorb ak and @k in the definition of the constant vectors ak and bk respectively but 
we leave them explicit on purpose. They play a very important role: for a given set of a\~ and bf~ 
we can tune the and (3k appropriately to ensure that the matrices Me(w) and M^ l (w) have 
unit determinant. Despite this, there is an ambiguity in the factorization (4.10) related to the 
normalization of the vectors and Nothing must depend on this choice of normalization. 
This will indeed be the case as will be apparent shortly. At this stage we just remark that the 
constant vectors at are the analog of the constant vectors of (3.16) in the BZ method. The 
ambiguity related to the factorization of rank one matrices in vectors is directly related to the 
ambiguity in the normalization of the vectors in the BZ method (cf. footnote 7). As is well 
known in the BZ method, nothing depends on the overall normalization of those vectors. 

In order to factorize A4e(vj) as in (4.1) we have to change from the constant spectral pa- 
rameter w to the space-time dependent parameter t through (2.20), which implies 

1 „JJ^ + i ), (4 .n, 



w-w k \t-t k 1 + iifc 

where the moving poles tk are determined by (2.20) evaluated at w k with the plus sign and 

= ~ h 2 lV (4-12) 

As a function of t, A4e(£, x) has a total of 2N poles: N poles at t = and N poles at t = — 1/tk- 
These have to be distributed among the factors A + and j4_ in (4.1). The analytic properties of 
the Riemann-Hilbert problem are such that the poles at t = —\jty~ come from A + (t) and those 
at t = t k from A_. One therefore makes the ansatze [19] 

N T 



fc=l 

N 



fc=i 



where the second equation arises in the factorization of M.^{w). These two equations introduce 
two new sets of vectors that we call c\~ and d k - 

The vectors a k , bk, Ck and dk are not all independent and determining their relation amounts 
to solving the Riemann-Hilbert problem. We first note that the constant matrices Ak and Bk 
are not independent since the two matrices A4e(il>) and M^{w) are inverses of each other. 
This determines the vectors bk from the up to scaling, a freedom that is reflected in the /3k 
in (4.9b). 

We can use the pole structure of M^t, x)M.^(t, x) to deduce some properties of the ak 
and bk- To start with, the absence of double poles at t = — 1/tk in the product implies that the 
vectors ak and bk are orthogonal: 

a\bk = for each k. (4.14) 
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From the absence of single poles at t 
the relations 



in the product M.^(t, x)A'i E 1 (t, x) one deduces 



a k a k a k A k 



A k b k (3 k b^, 



with the definitions 
A k 



1 + tt k 



A k 



M E (t,x) 



a k v k a k a[ 
l + tt k 



(4.15) 



(4.16) 



Equation (4.15) is satisfied if there exist j k such that 

alA k = ~i k v k fi k bl and A k b k = j k a k v k a k . 



(4.17) 



Then both sides of (4.15) are equal to a k i/ k f3 k r f k (a k b]^). We note that the j k defined uniquely 
by (4.17) depend on space-time. 

The next non-trivial step is to determine the vectors c k . With the ansatz (4.13a) for A + (t), 
the requirement that A + (t)Ai^ i 1 (t, x) have no poles at t = —l/t k gives the vector equation 

:, (4-18) 



c k T 



kl 



where T k i is an iV x N matrix with elements 



kl 



1L 

th 

t£n a k b i 



for 
for 



k = I 
k^l. 



(4.19) 



Solving equation (4.18) for the c k we obtain the matrix A+(t). A similar argument is used to 
construct A^^t) in (4.13b). One finds the equation for the vectors d k to be T k idi = a k . Solving 
this equation for the d k we can readily construct the monodromy matrix Me by taking the limit 
t — > oo in (4.1) and using that A^e(oo) = A + (oo) = 1. The result is 



A? 



M E = A+^oo) = 1+^ ftfctfc^r- 1 : 



kl a l 



(4.20) 



Finally we factorize Me = V^V^ and obtain the space-time fields. If needed, we can also 
construct explicitly the generating function Ve(£, x) from equation (4.3). 

At this stage it is instructive to investigate how the final answers (4.13a) (in conjunction 
with (4.18)) and (4.20) are insensitive to the ambiguity related to the rescaling of the vectors. 
Note that if we rescale the vectors a k to r k a k and b k to s k b k , then we must rescale a k and j3 k 
as r~^ 2 a k and s^ 2 f3 k respectively in order to preserve the products (4.10). This then means that 
the 7^ scale as rkSklk- It then immediately follows that the matrix T k i scales as r k siT k i. The 
inverse matrix (r -1 )^; naturally scales with the inverse factor 1 sf 1 and hence we see that the 
final answers (4.13a) (in conjunction with (4.18)) and (4.20) are insensitive to such rescalings. 
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4.3 Conformal factor 



For solitonic solutions of the previous subsection the conformal factor /e can also be obtained 
in a closed form. The final result is 



where /cbm is an arbitrary numerical constant. This constant needs to be chosen appropriately 
in order to ensure certain physical properties (say, asymptotic flatness) of the final space-time. 
A derivation of expression (4.21) is given in appendix A. 

5 Example: Kerr— NUT solution 

In this section we present a concrete implementation of the method of section 4 by constructing 
the Kerr-NUT metric. The construction illustrates all the steps of the BM solitonic method. 
More complicated examples are certainly doable; we leave such a line of investigation for the 
future. 

5.1 Construction of general 2-soliton solution 

The main difficulty in constructing the general multi-soliton solutions using the BM group the- 
oretic approach lies in finding meromorphic matrices A^e(w) that satisfy the coset constraints. 
The analog of this problem does not arise in the approach of BZ because they relax this con- 
straint and consequently have to renormalize the resulting matrices. This works well for SL(2, ]R) 
but already for SL(3,M) it gives spacetimes that do not represent black holes. A clever solu- 
tion of this problem in the BZ approach was found by Pomeransky [31]. In this respect the 
BZ method supplemented with the Pomeransky trick remains the most effective and powerful 
method for constructing solutions of vacuum four- and five-dimensional gravity. For a concise 
review and further references see [33]. The Pomeransky trick works well for SL(n,lR) but has 
no known analog for other coset models. As mentioned in the introduction, despite the initial 
complications, the promise of the BM method lies in its generality; it can be taken over to other 
coset models. 

It turns out that for SL(n, R)e monodromy matrices with a maximum of two poles, it is 
rather straightforward to take into account the coset constraints explicitly. When more poles 
are present one can perhaps set up a recursive algorithm for finding the appropriate meromorphic 
matrices. We have not attempted this and leave this line of investigation for the future. Here 
we present a discussion of a two soliton SL(2,M)e matrix. The most general such configuration 
describes the Kerr-NUT solution as we show now. 

We start with the general form of A4~e{w) (compare (4.9a)), 



/I = k B M ■ Y\_(t k u k ) ■ detr 



(4.21) 



k=l 



M E {w) = 1 + 



(5.1) 



w — c w + c 
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where 01 and a 2 are two-dimensional column vectors. The poles can be chosen in this way by 
a shift of axis, see (2.19). Given the constant 2x2 matrix a = (ai,a 2 ) and the 2x2 matrix 
£ = a T a we must choose 

2c f 62 \ 
a = , a = diagjai, a 2 ), (5.2) 



detM -£ 



11 



in order to satisfy the constraint det Me(w) = 1- The matrix M E {w) is symmetric and is of 
determinant one, hence it is in the coset SL(2, M)e/SO(2)e- For the parametrization of the 
inverse A^e(^) -1 we can choose b = a£ _1 e and (3 = — adet£ with 



■-(:-.')■ 

Due to the scaling freedom for the vectors, we can choose without any loss of generality 



(5.3) 



oi = U )' " 2= (1 • (5 - 4) 



Thus 



1 + Ci 2 C1 + C2 

2 
2 



C1 + C2 1 + C: 2 



(5.5) 



and 



2c(l + C|) 2 C (l + d 2 ) 

ai = jz — T2 «2 = -jz -T-TT2- 5.6) 

(1 - C1C2) 2 (1 - C1C2) 2 

For the bk vectors we have 
and for /3k 

/3 1 = -2c(l + C 2 2 ), /3 2 = 2c(l + Ci 2 )- (5-8) 

From the above expressions we see that afbi = = aTb2, as expected. Furthermore we have 
a^bi = — ajb2 = 1. With the above choices we find 



1 62 1+M2 

to - U \ 1 gl 2 t 1 (t 2 +t 2 1 ) 

\ 1 6l l+tlfe 



(5.9) 



From r we obtain the c& vectors by (4.18) and from there by looking at the limiting value of 
the AT 1 ^) at t = 00, cf. equation (4.20), we obtain M^ix). We find 



M E (x) = 1 + a(rTe- 1 0" 1 a T , where T = diag {ti, t 2 }. (5.10) 
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We also observe that 



TTe = {av) 1 — - — — — £, where (av) = diag {a\V\, a 2 v 2 }. (5-11) 

1 + t\t 2 



From Me(x) one can read off the physical fields. The conformal factor, given by (4.21), is 



2 t 1 Vit 2 V 2 
JE = fc BM" 



(t 2 -tl) 2 

5.2 Interpretation as Kerr— NUT metric 



& (i + q)(i + <3) ! 
6162 (I + M2) 2 



(5.12) 



For four-dimensional vacuum gravity with SL(2, M)e symmetry we can parametrise the mon- 
odromy Me as [3] 

*'-VfV--( A+ A A .7 (5.3) 

Here, ^ is dual to the metric function ^ by the duality relation 10 

A 2 

* 2 dil, = dip. (5.14) 

P 

The D = 4 metric is given by 

ds\ = -A (dt + ^dc))) 2 + A _1 (/| (dp 2 + (fc 2 ) + p 2 # 2 ). (5.15) 

(The form of the conformal factor is due to the change from Ehlers to Matzner-Misner variables 
to describe the physical space-time.) 

To write explicit expressions for the scalars it is convenient to introduce prolate spheroidal 
coordinates (u, v) 

z = uv, p = \J (u 2 — c 2 )(l — v 2 ), c < u < 00, — 1 < v < 1. (5.16) 
These coordinates allow us to write the pole trajectories t\ and t 2 as 

= (u-c)d + v) = (u + c)(l+v) 

yV -c 2 )(l -v 2 ) ^(u 2 -c 2 )(l-v 2 ) y ' 

The inverse relations are 

U u — c 1+v , 

-i = ^, ht 2 = - . 5.18 

t 2 u + c 1 — v 

In these new coordinates we have 

A = 2) \y 2(? ^ + & 2 + u2{l - ^ 2 -° 2{ - l+ tfx 1 + cl)] (5-i9) 



10 From the point of view of the Geroch group, this duality with its non-linear prefactor is at the heart of the 
infinite-dimensional symmetry and integrability of the system: ip is the first in an infinite set of so-called dual 
potentials on which the infinite symmetry acts. 
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and 

i> = 1 [2 OT (C 2 - Ci)(l " C1C2) - 2cMCi + C 2 )(l + C1C2)] , (5.20) 
where the common denominator of these expressions is 

D = v 2 c\Ci + C2) 2 + 2vc 2 (Cl ~ Ci 2 ) + n 2 (l - C1C2) 2 + c 2 (l + ( 2 )(l + Cl) + 2cu(l - CiCf)- 

(5.21) 

By applying the duality relation (5.14) one can write an expression for ip. n It is slightly more 
complicated 

* = O^V (5 - 22a) 

= u\l - C1C2) 2 - c 2 (l + d 2 )(l + Cl) " cV(Ci + C2) 2 , (5.22b) 
= -4c 3 Ci(l + Ci)(l + Cl) " 2c 2 (Ci + C 2 )(l - ClCDu + 2c(Ci - C 2 )(l - C1C2) V 
- 2c 3 (Ci - C2XI - CiC 2 ) 2 ^ + 2c(d - C2XI - C1C2) Vu 

+ 2c 3 (Ci + Ca)(l + 2Ci 2 + Ci 2 C 2 > 2 + 2c 2 (Ci + C2XI - CiCl)™ 2 . (5.22c) 

These expressions look somewhat cumbersome. To compare them with the corresponding ex- 
pressions for the Kerr-NUT metric in the standard Boyer-Lindquist coordinates, let us recall 
that the latter takes the form (see e.g. [4]) 

ds 2 = -1 (S - a 2 sin 2 0) dt 2 + ^dr 2 + UdO 2 + ^ ((£ + a X ) 2 sin 2 9 - x 2 H) # 2 

+ - (xB - a(E + ax) sin 2 0) dtd(j>, (5.23) 

with 

E = r 2 + (n + acos6>) 2 (5.24a) 

E = r 2 - 2mr - n 2 + a 2 (5.24b) 

X = asm 2 6 - 2n(l + cos6l). (5.24c) 

The Boyer-Lindquist coordinates (r, 0) are related to the prolate spheroidal coordinates (see e.g. 
appendix G of [34]) simply as, 

u = r — m, v = cos9. (5.25) 

We find that with the identifications 

Ci = , C2 = , c=ym z + n z — a z , (5.26) 

a + n c + m 



11 Alternatively, one could construct the generating function Ve(£) and apply the algebraic Kramer-Neugebauer 
transformation (relating the Ehlers and Matzner-Misner description) to obtain Vmm(£) that directly contains 

i> [3]. 
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the tt, t(f>, and (jxj> part of the metric matches with corresponding expressions obtained through 
(5.19), (5.20) and (5.22). Relations (5.26) can be inverted to read 

Cl + C2 Cl — C2 1 + ClC2 /_ n „N 

a = — m- — — , n = m- — — , m = c- — — . (5.27) 

I + C1C2 I + C1C2 I-C1C2 V ' 

When £1 = ^2 we obtain the Kerr solution and the corresponding expressions match those of 
[25] and [19] (when certain minor typos and misprints are fixed in those references). When 
(1 = (2 = we obtain the Schwarzschild solution as in [25, 32]. 

The conformal factor can also be easily computed using the formula (4.21). We find 

2 , u 2 — m 2 — n 2 + a 2 v 2 

Je = ~ &BM77 — q— ■ — 2v 2 2~~2V (5.28) 

4(m/ + n z )(u z — c z v z ) 

Choosing the constant /cbm to be — 4(m 2 + n 2 ) so that /e — > 1 as r — > 00, we have 

= U »-m»-n'H-V 

(u z — C Z V Z ) 

which in the MM coset allows us to match directly with the metric (5.23) 

A-/! = (m + 1 2 + ( 2 n 2 +mj)2 . (5.30) 

For completeness, we record the form of the (Ehlers) monodromy matrix A^e(^) expressed 
in terms of the physical quantities 



(m + w) 2 + (n + a) 2 2(am-nw) 

Me(w) = 5- ' ' . \, , (5.31 

w z — c z \ 2 (am — nw) (w — m) + (a — n) 1 



with c as in (5.26). 



5.3 BZ Ehlers Construction 

The Kerr-NUT metric has been reconstructed by several authors in the context of the standard 
BZ method, see e.g. [4, 33]. In this section we revisit this computation and perform it using the 
Ehlers generating function 1 3/e(A). To the best of our knowledge this has not been presented 
before. 

The Killing part of the metric of flat space translates into Me being the identity matrix. 
This implies that the seed \^e(A,x) is also identity. Now we add solitons at A = \i\ with BZ 
vectors m.Q = {A\,A-i) and at A = \i<i with BZ vectors = (Bi,B2). The Tbz matrix is 
found to be 

/ A\+A\ A 1 B 1 +A 2 B 2 \ 

and we need to rescale the resulting Me by — ^Jt 2 [2, 4] to ensure the determinant condition. 
With the choice of parameters w\ = — c and W2 = +c for the pole locations /zi and fi2 in (3.15), 
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and (Ai,A2) = (— Ci> 1) an d (-Bi, -B2) = (1) — C2) for the BZ vectors, the final expression for the 
rescaled Me matches precisely with the Me obtained from the BM method in section 5. The 
conformal factor obtained via equation (3.22) also matches with the one obtained in section 5, 
with the constant k^z in (3.22) being k^z = 4(m + c) 2 . It is intriguing to note that up to overall 
normalizations the BZ vectors are precisely the vectors b\ and 62 of the BM construction above. 

6 Conclusions 

In this paper we studied the integrability of two-dimensional gravity-matter systems with matter 
related to a symmetric space Ge/Ke- The integrability of these set-ups can be exhibited either 
using the Belinski-Zakharov (BZ) linear system or through the Breitenlohner— Maison (BM) 
linear system. As emphasized, the approach of Breitenlohner-Maison makes the group structure 
of the Geroch group manifest. We analysed the relation between the BZ and the BM linear 
systems and presented explicit relations between the generating functions appearing in these 
two linear systems. 

An embedding of the Belinski-Zakharov solution generating technique in the Geroch group 
was also studied in general. We pointed out that it is impractical to find a satisfactory general 
embedding of the full BZ solution generating technique in the Geroch group. Relation (3.10) 
provides in principle the link between the generating functions. However, one must keep in 
mind that the left hand side of (3.10) must be a 'physical' generating function ^e(A). Here, 
by a 'physical' generating function we mean a generating function that gives a representative of 
the coset Ge/Ke upon taking the limit A — > 0. This does not happen automatically in the BZ 
technique where the dressing matrix xW is considered in a more general context. 

On the other hand, following the unpublished work of Breitenlohner and Maison [19], we 
exhibited a novel solution generating method where the group theoretical interpretation is clear 
from the beginning to the end. In our approach we solve the requisite Riemann-Hilbert problem 
algebraically. Since only algebraic manipulations are involved, our technique is akin to the BZ 
technique. As a novel example, we constructed the Kerr-NUT solution in this approach. 

Our main interest in performing the analysis of integrability in these gravity-matter systems 
is to make the means for constructing new solutions available in situations where the standard 
inverse scattering method of BZ is not applicable. This is typically the case for extended 
supergravity theories that have a string theory origin. For minimal D = 5 supergravity with 
exceptional Ehlers symmetry Ge = G2(2) this problem was pointed out in [15] and also arises 
for the STU model [35] with Ge = SO(4,4) or maximal supergravity with Ge = Es(8) [36]. The 
method explained in section 4 is still applicable in those cases as long as one finds a way to 
parametrise the Ge valued monodromy A^e(^) in a way similar to (5.1) and (5.2). There are 
several ways in which this could be achieved: (i) One could use global elements k S Ke C Ge C 
(Geroch group) to rotate the vectors in A4(w) in (5.1) into canonical positions. In the example 
of section 5 this would correspond to setting (1 = £2- The solution for canonical vectors can 
then be generalised by applying conventional Ke charging transformations. These charging 
transformations use only a very small subset of the full power of the Geroch group, (ii) One 
could embed Ge in GL(n,R) for n large enough and then solve the constraints on the vectors 
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for the embedding explicitly. We expect a combination of these two techniques to be the most 
promising line of attack. Alternatively to (i) and (ii), one could perhaps hope to develop 
some general algorithms to find the appropriate monodromy matrices by combining ideas from 
uniqueness proofs for black holes, see e.g. [37], and the fact that the monodromy matrices are 
closely related to the behavior of solutions on the z-axis, see e.g. section 4 of [3]. In future work 
we plan to explore these sets of ideas. See also [30] for a slightly different but related viewpoint 
on this problem. It will be very interesting to see if the algebraic Riemann-Hilbert factorization 
approach can be used to construct new black hole solutions generalizing [20, 21] and new fuzzball 
solutions generalizing [22, 23]. 

Another interesting aspect of the integrable structure in two-dimensional models is their pos- 
sible relation to the recently studied infinite-dimensional symmetries of string and M-theory [38, 
39]. These symmetry groups are extensions of the Geroch group and are conjectured to be 
symmetries of the unreduced theory. First steps in investigating this relation were undertaken 
in [40] in the restricted case of polarized Gowdy space-times. The relation (2.18) between a spec- 
tral parameter and a space-time coordinate suggests that a mapping between Lie algebraic and 
geometric data might be possible. If taken seriously, this approach would allow the treatment 
of space-time as a concept that fully emerges from symmetry considerations. 
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A Computation of the BM conformal factor 

In this appendix we present a derivation of equation (4.21) for meromorphic monodromy matrices 
following [19]. Using the light cone coordinates (2.13) with the property *2<9± — ±id± we can 
write the differential equations for the conformal factor f E as 



Next we wish to write Tr(i-E,±-PE,±) in terms of the matrix A + (t) introduced in (4.3) in section 
4. To this end we evaluate the residue of the poles at t = rti in the Lax equation (2.17). For 
evaluating the residue on the l.h.s. of (2.17) we use the relation 




(A.l) 



d±v E (t,x) = d±v E (t,x)\ t + (d ± t)v E (t,x) 



(A.2) 



where V E (t, x) 



dV E (t,x) 
at 



. These relations together with (2.18) give 



±id±lnpV E {±i) = P E ,±V E {±i). 



(A.3) 
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Now replacing (4.3) in (A. 3) we obtain an expression for Pe,± in terms of A+(±i) and A + (±i). 
Substituting that expression in (A.l) we obtain 

9±ln/ E = ~(d ± kip)Tr (^(iiJA^ii)) 2 . (A.4) 

Using the explicit form of A + (t) from equation (4.13a) together with (4.18) and the identity 
a^bk = T k i(t k - ti) we get, 

A7 1 (t)A+(t) = -b 1 T" 1 1 a T . (A.5) 

In writing this equation we have used a convenient matrix notation, where T is a diagonal matrix 
with entries if.. Differentiating for k ^ I the identity (cf. (4.19)) 

a T k bi = T k i{t k -ti) : (A.6) 

with respect to the light cone coordinates we obtain the components with k ^ I of the equation 

d±T = -(5± hip)-^- [V T iTT T iYT + TFT] — ?— . (A.7) 

Looking at the definition of the diagonal components in (4.19), we note that we need d±^k 
in order to obtain the corresponding expression for the diagonal components of V. We hence 
differentiate the relation A^bk = 7k®kVkak from (4.17). An important intermediate result for 
this is that Akbk is constant. From this it is easy to deduce that (A.7) holds as well for the 
diagonal components. 

Using these formulas, we first substitute (A.5) into (A.4) and then manipulate the new r.h.s. 
to bring out terms that are total derivatives using (A.7). As a result, we can rewrite equation 
(A.4) in the form 

3 ± ln/ E = ilr (r-^r) + ^Tr ((Tv)~ 1 d±{Tv)) , (A.8) 

where {Tv) is the diagonal matrix with entries t^v^. Equation (A.8) can now be readily inte- 
grated to give the final result (4.21) 

N 

/i=fc B M-n^)- detr ' (a - 9) 
fc=i 

with /cbm an integration constant. (More generally, the conformal factor is related to a cocycle 
calculation in the affine group [3].) 



B BZ Matzner-Misner 

The Belinski-Zakharov (BZ) approach is a well established solution generating technique for 
vacuum gravity. The method is applicable in any dimension, though only in four and five 
dimensions can the generated solutions be asymptotically flat [33]. In this appendix we focus 
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on D = 4. We assume that the space-time admits two commuting Killing vectors, one spacelike 
(angular) and one timelike. In this case the four-dimensional metric admits the following form 
in the Weyl canonical coordinates (cf. (5.15)) 

ds 2 = e 2u (dp 2 + dz 2 ) + g ab dx a dx b , (B.l) 

where the indices a, b run over the Killing coordinates cp and t. The vacuum Einstein equations 
can be used to choose without any loss of generality the coordinates such that [34] 

detg = -p 2 . (B.2) 

For this class of metrics, the Einstein equations divide in two groups: one for the Killing 
part g of the metric 

8 p U + d z V = 0, with U = p(d p g)g~ l and V = p(d z g)g l (B.3) 

and the second group for the conformal factor v. From the discussion in the main text of the 
paper it is clear that equation (B.3) are the equations for a GL(2,K) integrable sigma model. 
The equations for the conformal factor v read 

d p v = -^ + ^-Tx{U 2 -V 2 ), d z v = ^Tr(UV). (B.4) 

2p 8p Ap 

Note that the equations for g do not contain the function u. Once one obtains a solution of 
(B.3), v can be obtained by a line integral. 

The Belinski-Zakharov spectral equations for the above GL(2,R) model (B.3) are [2] 12 

*. = £^., *. = f±£., 

where A is the spectral parameter and D\ and Di are the two commuting differential operators 
of (3. 3). 13 The generating matrix ^(A, p, z) is such that in the limit A — > it gives the volumeful 
metric g. Using the above linear system (B.5) one can construct an infinite class of new solutions 
by dressing seed solutions. This procedure has been reviewed at several places, see e.g. [4, 33]. 
The matrix g can also be written in terms of the SL(2,M)mm unimodular matrix Mum 

g = pM MM , (B.6) 



where the subscripts MM stand for Matzner-Misner. In terms of the coset variables it takes the 

'mm 1 



form Mmm = VmuV^mm, where Vmm is the Matzner-Misner coset representative of the quotient 



12 Note that there is neither an Ehlers nor a Matzner-Misner subscript on this generating function since it 
agrees with neither. It is, however, closely related to the Matzner-Misner version as will become clear in the 
sequel. 

13 In 'light-cone' coordinates the spectral equations read 

D ± 9= ±if '[ a ± 9 ) g ~ 1 9, with D± =d ± --^d x . 
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SL(2, R)mm/S0(1, 1)mm an d r] is the S0(1,1)mm invariant metric rj = diagjl, — 1}. Following 
[3] we take the parameterization for Mmm to be 



M MM = V£ m t)Vmm = ( A _/^ _% J with r, = I Q ^ J , (B.7) 

where we emphasize that in our conventions the second index denotes the time component. We 
use this (non-standard) convention to facilitate comparison with [3]. 
Equations (B.3) and (B.4) can now be brought to the form 

d m (pM M ^ m M MM ) = 0, (B.8) 

and 

= | (Tr (M-^dpMuu) 2 ~ Tr (M m ^ 2 M M m) 2 ) , (B.9) 



= ^Tr {M^ M d p M MM M^ M d z M MM ) , (B.10) 
where we have used Dp = and £(p, z) is defined as £ = e v . Defining as before 

flviM,m = r (9 m 14lM^MM + (^iti^MM^mm) ) J (B-H) 

we can rewrite the above equations in the coordinates x^ of (2.13) as 

±iC l d±ti = ^Tr (P M m,±Pmm,±) , (B.12a) 
AnO^MMHO- (B.12b) 

Equations (B.12a) and (B.12b) are formally identical to equations (2.14a) and (2.14b). The BZ 
Lax pair can be written for these equations [15] as well. 

If one wants to relate the BZ-generating function ^(A) to the group-theoretic BM-generating 
function Vmm(£) additional care has to be taken because of the factor of p in (B.6). A convenient 
choice is 

*(A, x) = v / 2pl^V^ M (x)riVMM(t, x). (B.13) 

(Note that this differs from what was given in [3]; their choice does not map the linear systems 
into each other away from t = A = 0.) 
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